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Abstract 



We report on a generalized Wulff construction that allows for the calcu- 
lation of the shape of two-dimensional materials with orientational order but 
no positional order. We demonstrate that for sufficiently large domain radii, 
the shape necessarily develops mathematical singularities, similar to those re- 
cently observed in Langmuir monolayers. The physical origin of the cusps 
is shown to be related to the softness of the material and is fundamentally 

diffferent from that of the sharp angles seen in the shape of hard crystals. 
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The singularities in the shape of crystalline materials have long been understood to be 
a macroscopic expression of the positional order of crystals at the atomic level. The angles 
between the faces of a crystal are, for instance, structural invariants dependent only on a 



construction that allows one to determine the equilibrium shape of a crystal, provided one 
knows the values of the surface energies for the various Miller indices. If, for every index, 
the corresponding surface is placed a distance from a fixed point proportional to the surface 
energy, then the inner envelope of the planes constitutes the minimum energy crystal shape. 

Surprisingly, sharp edges in the shapes of samples are not just encountered for hard, 
crystalline materials. Using the Wulff construction, Herring argued that liquid crystals 
without positional order, but with orientational order could also display sharp ables, or cusps 
0. The cusps are no longer material invariants. Cusped domain shapes are, indeed, well- 
documented for three dimensional liquid crystals ||||. More recently, studies of the shapes of 
two dimensional materials without positional order have also revealed cusps and sharp angles 
These studies have also revealed that the internal structure of such "soft" materials is 
significantly deformed and dependent on the domain shape, while Herring assumed a rigid 
internal structure. Thus, the Wulff construction is not manifestly valid in this case. 

In this paper we present a formalism which allows for the construction of domain shapes of 
soft two dimensional materials when the internal strucure is described by a two dimensional 
XY model (e. g. hexatic or nematic liquid crystals). We will demonstrate that cusps are, 
indeed, not just possible, but that they ought to be a generic feature of such domain shapes. 
Moreover, the cusp angle provides important information on the elastic moduli and surface 
energy of the material. 

We will model the internal structure of a domain by a unit vector c = (cos 0(f), sin 0(f)). 
The associated XY model free energy is 



Here, k is the stiffness (or Frank constant) of the order parameter field and a is the 
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certain set of integers (Miller indices In a classic paper ||, Wulff developed a geometrical 




(1) 



anisotropic surface energy, which depends on the relative angle, G — 9, between the or- 
der parameter and the outward normal h = (cos 9, sin 9) to the domain boundary. Note that 
a(0-2vr) =<r(0). 

To find the optimal domain shape we must minimize H [0(x,|/)] with respect to Q(x,y) 
and the domain shape, while keeping the domain area, A, fixed. We start with the limiting 
case k = oo when the texture is rigid so that = and only the shape needs to be varied. 

K = OO 



The domain shape is a closed planar curve. We will express the co-ordinates (x, y) of a 
point on the curve in terms of the angle 9 of n and the minimum distance R{9) between the 
tangent line through (x,y) and the origin. In terms of R{9): 



x{9) = R{9) cos(#) - < ^p- sin(0) 

d9 

y (9) = R(9)sm(9) + ^p-cos(9). 
The variational equation determining the domain shape is, then 
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(3) 



The first term above is just the surface energy, with (R{9) + d 2 R(9)/d9 2 ) d9 a line element 
along the domain boundary. The second terms is — XA, with A a Lagrange multiplier. The 
variational equation reduces to 



R{9) + 



d 2 R(9) 
d9 2 
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A 



a{9) + 



d 2 a{9) 
d9 2 



(4) 



with R(9) = <j(9)/\ as the solution, modulo an overall translation of the domain. The 
domain shape R(9) =a{9)/\ corresponds precisely to the Wulff construction. The somewhat 
unusual parameterization of two dimensional curves thus allows for a straightfoward and 
analytical determination of the domain shape. 

The results above have been derived previously J7|. However, it has not, to our knowledge, 
been noted that variational derivatives with respect to R{9) can be carried out when the 



anisotropic surface energy of an element of boundary depends on the location of the element 
as well as its orientation. One simply replaces a (6) by o~(x, y, 9) and parameterizes x and y 
as in Eqs. (0) . In the remainder of this Letter we describe the consequences of this strategy 
as applied to a two-dimensional domain. 

k finite 



We now allow the internal structure of the domain to respond to changes in the domain 
shape. Minimizing H [0(x, y)\ with respect to the angle Q(x,y) leads to the requirement 



V 2 e(x,y) = 



de(x,y) 



dn 



- o-'{e-Q{x,y))\ houadaxy 







(5a) 
(5b) 



boundary 

In Eq. (5b) the d/dn derivative is along the outward normal. With no loss of generality we 
can write the solution of Eq. (5a) as 



©(z> y) = ~ (f(x + iy) - f(x - iy)) , 
i 



(6) 



with f(z) an arbitrary analytic function. It is convenient to rescale x and y by the mean 
domain radius Rq- The domain boundary in the complex plane of a nearly circular domain 
is, then, the unit circle z = e td . Eq. (5b) then reduces to 



^ (zf'(z) - l/zf'(l/z)) = io< {-% log^ + i (/(*) - /(!/*))) 



(7) 



with z on the unit circle. 

If we now go through the same free energy minimization as for k = oo, we find instead 
ofEq. §) 

d 2 R(8) 



R(6) + 



dO 2 



(8) 



with o"o an isotropic surface tension and 
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-[f(z)+f(l/z)) 



+ {zf(z)Y+(-f(l/z)) + I dw 



f'M 2 - -JW 



(9) 



Equations (^) and the latter relation holding if the domain is not too deformed from a 
perfect circle, are our key results. If we know the function f(z) on the unit circle — through 
Eq. (0) — then we can reconstruct the domain shape with the use of Eqs. (§) and Note 
from the discussion following Eq. (^) that we can interpret the solution R(9) of Eq. (ffi) as 
being proportional to the effective surface tension, i. e. a surface tension that incorporates 
the interior softness of the domain. 

Finding the complex function f(z) looks like an intractable problem, as Eq. (|7|) is highly 
nonlinear. We will consider some physically relevant forms for a ((f)) to show how f(z) and 
the domain shape may be found. 

cr((f>) = <j + a n cos n0 

The case n = 1 corresponds to an anisotropic surface energy proportional to c • n. The 
case n = 2 corresponds to the anisotropy energy of a two dimensional nematic. 
It can be verified, by direct substitution, that 

/(*) = ^ log (1 (10) 

is a solution of Eq. (|7|), with 

a n R S = l " a " R " /K (11) 

1 + Y 1 + (na n R / k) 

When n = 1 the texture corresponds to the two-dimensional version of a "virtual boojum" 
HI! — a singularity in the texture that lies outside of the domain. Imaging of domains of 
monolayer and near-monolayer films by Brewster angle microscopy reveals strong evidence 
for the existence of this structure. ||10|-[T2|. The distance of the boojum from the center of 
the domain, Rb, is given by 



Ss ^Wwf. (12) 

a\RjQj k 

If the parameter V = n/aiRo is small compared to one, then Rb approaches Rq, while in 
the limit r = oo the boojum retreats to infinity. The texture is, clearly, highly deformed 
when r << 1. 



The domain shape is found by substituting from Eq. (|10D into Eq. @. Surprisingly, 
F{0) = 0, so R{9) = -Ro- The domain shape is thus a perfect circle for n — 1. For n > 1 the 
exact solution for the texture is equivalent to the texture produced by n singularities lying 
outside of the domain. These singularities are no longer boojums in the accepted sense of 
the term |||9| . The order parameter angle advances by 47r/n as one traces a path encircling 
one of the singularities; an advance of 4n characterizes the boojum. The singularities for 
n > 2 are "fractionally charged" in the topological sense. However, the fractional nature of 
these singularities appears to have no consequence for the textural structure of interest, as 
the singularities lie outside of the domain. Now T{9) is no longer equal to zero. The domain 
shape is deformed from perfect circularity. 

O"(0) = CT + a l cos + a 2 COS 20 



This is the anisotropic boundary energy believed to be relevant to domains of liquid- 
condensed phase in Langmuir monolayers. It represents the lowest three terms in a system- 
atic Fourier expansion of the anisotropic surface energy. For 02 « ax, the function f(z) is 
given by 

/(,)= log (!-«,) -3 ° 2 ;~ a /'fW^> \£z!L] dt + o(( a -l) 2 ), (13) 

a\ 1 — a 1 — az Jo 11 — aztl \\ai/ J 

as can be checked by direct substitution. The small admixture of 02 cos 20 in cr(0) has 

dramatic effects on the shape. Using Eq. (|I3D in Eqs. (||) and (pi), we find that 



For R < 0.356<7oft/aia2, the domain has a smooth, nearly circular, shape. 

For Ro > 0.356o"ofi;/aia2, the domain has a cusp singularity in its shape. For i? - ► 00, 
the cusp angle Aip — the difference between the inner angle of the cusp and 7r — obeys 

A^(R ) = 5.48 K/aiRo + O (l/R 2 Q ) . (14) 

Note that the asymptotic variation in the excluded angle is independent of the pa- 
rameter 02- The O2cos20 term is thus a singular perturbation. The domain remains 
nearly circular outside of the immediate vicinity of the cusp. 
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Large domains will, thus, always have a cusped boundary. Figure [T] graphs the evolution 
of the excluded angle of the cusp as a function of the radius of a domain in which surface 
and bulk energy parameters have values that are consistent with the expansion described 
above. Figure ^ depicts a cusped domain as generated by the generalized Wulff construction 
for a specific domain radius. The inner angle of the cusp, ?/v, is indicated in the Figure. 

The measurement of cusp angle can be achieved by visual inspection of the domain. In 
light of the results reported above, a plot of versus Rq should yield the parameter ratio 
n/ai, and also the combination aon/a^. A measurement of the cusp angle as a function 



of domain radius has, in fact, been carried out Hl3[ . The results will be reported elsewhere. 

In summary, we have found that deformable domains of materials with an ATY-like order 
parameter may have shape singularities, even though their texture is perfectly analytic in 
the domain of definition. It is important to note that the material softness plays a key role. 
If we set the stiffness, k, to infinity, then the boundary does not have a cusp. The physical 
origin of our cusp is thus strikingly different from that of hard crystalline materials. It would 
be interesting to know whether the results reported here extend to three dimensions. 
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Theorique at Saclay, where some of this work was carried out. 



7 



REFERENCES 

[1] See, for example, C. Kittel Introduction to Solid State Physics, 6th ed. (John 
Wiley and Sons, New York, 1886). 

[2] G. Wulff, Z. Krist 34, 449 (1901). 

[3] C. Herring in Structure and Properties of Solid Surfaces, R. Gomer and C. Smith, 
eds. (University of Chicago Press, Chicago, 1953). 

[4] G. Friedel and F. Grandjean, Bull Soc. Fr. Mineral 33, 409 (1910); J. B. Fournier and 
G. Durand, J. Phys. II Fr. 1, 845 (1991). 

[5] P. Oswald and F. Melo, J. de Physique 50, 3527 (1989). 

[6] H. Bercegol, F. Gallet, D. Langevin and J. Meunier, J. Phys. Fr. 50, 2277 (1989). 

[7] W. K. Burton, N. Cabrera and F. C. Frank, Philos. Trans. Roy. Soc. 243, 299 (1951). 

[8] S. A. Langer and J. P. Sethna, Phys. Rev. A 34 5035 (1986). For the original discussion 
of boojums, see ref. 0. 

[9] N. D. Mermin in Quantum Fluids and Solids, S. B. Trickey, E. Adams and J. Duffy, 
eds. (Plenum, New York, 1977). 

[10] D. Honig, G. A. Overbeck, and D. Mobius, Adv. Mater. 4, 419 (1992). 

[11] J. Meunier, unpublished. 

[12] See also the polarized reflection micrograph of a freely suspended liquid crystal film by 
N. A. Clark, D. H. Van Winkle and C. Muzny in ref. §. 

[13] D. Schwartz and C. Knobler, unpublished. 



8 



FIGURES 

FIG. 1. The excluded angle, Aip, of the cusped domain, defined as the difference beween ip- m , 
the inner angle of the cusp (shown in Figure 2), and it. The plot is of the excluded angle versus 
domain radius when two of the parameters controlling surface and bulk energy satisfy 02/0" = 0.05. 
The domain radius is in units of n/a\. 

FIG. 2. the shape of a domain when the parameters controlling surface and bulk energy are 
such that the domain has a cusp. Note that the domain is nearly circular, except in the immediate 
vicinity of the cusp. Barely visible in the Figure is the "swallow tail" that appears as an appendage, 
attached to the domain at the cusp, when one implements the Wulff construction. This appendage 
is amputated to generate the true domain shape. The virtual boojum lies outside of the domain 
and near the cusp. In the Figure, 02/0" = 0.25 and Rqcli/k = 5. 
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